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FRIS PROJECT OVERVIEW 

The Forest Resource Information System Project (FRIS) Is a coopera- 
tive effort between the National Aeronautics and Space Administration 
(NASA) and St. Regis Paper Co. (STR). Purdue University's Laboratory 
for Applications of Remote Sensing (LARS), under contract to NASA, will 
supply technical support to the project. 

FRIS is an Application Pilot Test (APT) Project funded by NASA. 

The project is interdisciplinary in nature involving expertise from both 
the public and private sectors. FRIS also represents the first APT to 
involve a large broad base forest industry (STR) in a cooperative with 
the government and the academic communities. 

Purpose 

The goal of FRIS is to demonstrate the feasibility of using computer 
aided analysis techniques applied of Landsat Multispectral Scanner Data 
to broaden and improve the existing STR forest data base, thereby 
creating the foundation of a dynamic information system. The successful 
demonstration of this technology during the first half of the project 
will lead to the establishment by STR of an independently controlled 
operational forest resource information system in which Landsat data is 
expected to make a significant contribution. FRIS can be viewed by the 
user community as a model of NASA's involvement in practical application 
and effective use of space technology. Additionally, FRIS will serve to 
demonstrate the capability of Landsat MSS data and machine-assisted 
analysis technology to private industry by: 
o Determining economic potentials, 
o Providing visibility and documentation, and 
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o The ability to provide timely Information 
and thus serve management needs. 

The ultimate long tern success fulness of FRIS can be measured through 
future development of remote sensing technology within the forest 
products Industry. 

Scope 

FRIS Is funded as a modular or Phase project with an anticipated 
duration of three years. The original project concepts were developed 
in 1973, and a formal project plan was submitted to NASA by STR In 
1976. The project officially began In October 1977 after the signing 
of a cooperative agreement between NASA and STR i and after the compel - 
tion of contractual arrangements with Purdue University. 

Organization 

The organization of FRIS is depicted in the chart that follows. 
Since FRIS Is a cooperative Involving three independent agencies, a 
steering committee consisting of a project manager from each Institution 
was formed to provide for overall guidance and coordination. Opera- 
tionally, both STR and LARS have project managers and project staff to 
Insure for the timely completion of activities within the project. 

The NASA technical coordinator monitors project activities and provides 
a liaison between the STR and LARS staffs. The solid lines on the 
chart indicate the flow of management responsibility. The dash lines 
reflect the technical and scientific Inter-changes between operating 
units. 


FRXS Organisation 


Steering Commit tea 

ASVT Project Manager 
NASA Technical Monitor 
FR1S Project Manager 

f 1 1 

Resource and Technology ___ NASA LARS/ 

Department/ STR . Johnson Spacecraft Center Purdue University 


I — Computer Systems - - -- -- -- -- - • Systems Design 


L. Cartographic Systems 


Mapping Unit 


U Forest Sampling Systems 


Classification Unit A 


L Cost Analysis 


• Cost Unit 


1 


1.0 INTRODUCTION 

The material appearing In this report Is a reflection of the FRIS 
Project Staff activities for the period 1 October 1979 to 31 December 
1979. This time frame encompasses the third quarterly reporting period 
for Phase III of the Forest Resource Information System (FRIS) Appli- 
cations Pilot Test (APT). Phase III of FRIS Is directed at meeting the 
overall Project goal: 

To document and transfer remote sensing technology developed 
throughout the project that will provide St. Regis with an 
Independent operational system, having Landsat data as a 
significant and viable contributor. 

Significant strides we< • made toward reaching this Project goal during 
this reporting period. Noteworthy activities that occurred during this 
quarter are summarized below: 

o A terminal connecting LARS to St. Regis National Computer Center 
was installed. 

o LARS staff received training on use of the St. Regis terminal, 
o The System Design Committee defined functional specifications 
for the FRIS configuration. 

o Preliminary versions of new software was prepared for testing, 
o Results of the applications test initiated the previous quarter 
were presented to St. Regis. 

o A technique designed to predict age of planted pine stands was 
tested for the Knabb Application Test Area. 

The remaining sections of this report describe these activities in 


mort complete detail. Appendix B contains update timeline charts for 
all tasks. 
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2.0 TASK AREA ACTIVITIES & 

2.1 Technology Transfer Task 

A key element In the System Transfer Phase of FRIS Is the develop- 
ment of a knowledgable core of St. Regis analysts. Formal short 
courses, workshops and text material are Items which help support this 
development. However useful these exercised are, they are no sub- 
stitute for actual analysis experience. 

The “hands-on" yperlence requirement for analyst training was 
the motivation for supporting a remote terminal link between the 
Southern Tlmberlands Divisional office at Jacksonville, FI and LARS. 
Details regarding this remote terminal configuration and operation can 
be found in previous FRIS Quarterly Reports. Previous task activities 
were directly related to supporting this remote terminal link. There- 
fore, preliminary training and documentation supports the Interactive 
LARSYS analysis capability of the terminal. Operationally, the FRIS 
Image Processing will be conducted in a batch mode. This chanqe will 
necessitate modification of the training procedures and documentation 
to reflect the St. Regis National Computer Center (NCC) operating en- 
vironment. 

In order to accomodate the smooth transition of the technology 
associated with the change In computer operating environments a decision 
was made to connect LARS to NCC via a remote terminal link. This link 
would allow the Technology Transfer staff to develop supplemental 
training and documentation necessary to execute LARSYS functions through 
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ROSCOE, St. Regis* Interactive echtlng software. 

During this quarter, the equipment for Project staff terminal 
connection to Dallas was ordered, received and Installed. It was pre- 
viously determined that this terminal conncectlon would consist of an 
IBM 3275 CRT terminal for edtttlng and job entry (with the necessary 
phone line and modem) to be provided by St. Regis, and a Data 100 RJE 
terminal for receiving printer output (with the equipment required to 
provide dial-up access to the St. Regis computer), which would be 
supplied by Purdue/LARS. In October, the following equipment was 
ordered by Purdue/LARS to allow us to connect the Data 100 terminal 
at LARS/Flexlab 2 to the St. Regis NCC: a Racal/Mllgo MPS48 dial 
modem, a T-bar digital switch, two cables and a business telephone 
with a manual DAA. 

Also during October, the IBM 3275 CRT terminal arrived. Its 
telephone line was Installed In early November and the modem arrived 
November 12. By November 16 the terminal, modem and phone line were 
connected together and working. For greater security, this terminal 
has been placed In the Purdue/LARS computer room and the LARS computer 
operators have been given a list of project personnel authorized to use 
the terminal . 

The telephone for the Data 100 terminal connection was Installed 
In mid-November, the dial modem arrived November 29 and was Installed 
December 7. After a few more t hinges to the modem and the removal of 
a 4-wtre modification In the Data 10Q load deck, the Data 100 terminal 
was able to receive printer output from the St. Regis NCC. 


s 


Training Activities 

On Monday and Tuesday* October 8 and 9* A1 Beecher Mas at LARS to 
review hU progress on the Installation of the LARS Data Reformatting 
program on the St. Regis NCC computer. He discussed the LARSYS 
support routines with Bill Shelley, the reformatting routines with 
Chuck Smith and the geometric correction programs with Cathy Kozlowskl. 
No major problem are anticipated. 

A ROSCOE training session scheduled for November 27 was abbreviated 
to an Informal demonstration by A1 Beecher of sign on and edlttlng 
procedures, due to the delay In the Installation of the Data 100 con- 
ruction. 

The ROSCOE training sesrlon was rescheduled for December 19 at 
LARS. Bob Lynn presented the material to ten project personnel at LARS 
and assisted them In hands-on sessions at the IBN 3275 CRT terminal. 

2.2 System Transfer Tasks 
2.2.1 System Design Committee 

During the third quarter of Phase III the System Design Committee 
met several times to formalize FRIS specifications. The list of 
functional specifications that were developed by the committee appear 
In Table 2.2.1 . Three vendors of data base systems were asked to 
demonstrate their systems capabilities and bid on the system Installa- 
tion In Jacksonville, FL. 

Demonstration materials were prepared by FRIS staff. Each vendor 




Table 2.2.1 Functional Specifications for Evaluation of FRIS System 
Design Alternatives 

I. Graphic Data Capability 

A. Input 

B. Analysis 

C. Update 

D. G'tput 

II. Tabular Data Capability 

A. Input 

B. Analysis 

C. Update 

D. Output 

III. Image Data Capability 

A. Conversion from vector to grid 

B. Conversion from grid to vector 

IV. Other 

A. Hardware 

1. Configuration 

2. Deliverablllty 

3. Support 

4. Data Communications 

B. Software 

1. Availability/cost of source 

2. Support 

3. Transportability 

C. Implementation 

1. Cost 

2. Time 

D. Vendor Profile 

1 . Customer base 

2. Customer service 

3. Expertise In forest based applications 

4. Vendor stability 

V. Overall Cost 





Table 2.2.2 FRIS data base manipulation requirements 

1. Produce a plot of the digitized data, containing the AU (Adminis- 
trative Unit} and OA (Operating Area) boundaries for all four of 
the AU's. 

2. The fourth file of the tape contains some extraneous points, pro- 
duce a clean plot demonstrating the edlttlng capabilities. 

3. Convert the Landsat classification data from grid to vector format. 

4. Produce a plot of each layer of information 

a. AU boundaries 

b. OA boundaries 

c. Landsat classification 

5. Associate attribute data with each layer of information 

a. for the AU boundaries layer, the attributes would consist of 
the AU numbers (264, 267, 268, and 271). 

b. for the OA boundaries layer, the attributes of Interest would 
be the OA numbers, the forest type, and the age of the stand 
(this information may be found on the sheets describing each 
Individual AU). 

c. for the classification data, this would be the names of the 
classes taken from the classification results tape. 

6. Produce an overlay of the three layers of information. 

7. Graphically represent where the Landsat classification and the map 

are in disagreement for a cover type. What we have in mind is a 
map depicting areas that would satisfy such Boolean combinations 
as: NONSTOCK (from the Landsat classification) .AND. .NOT. (forest 

types 9 .OR. 92 (from the map)). 

8. It would also be desirous to have maps of areas based upon the 
attributes of the OA's (e.g., Forest types 2, 11, and 21 which are 
greater than 15 years old). 

9. demonstrate the capability to apply transformations to the vector 
data sets (e.g., for rotation and scale). 


received the following data: 

1. Map of All's (Administrative Units) 264,267,268, and 271. 

2. Documentation of map contents. 

3. Tape containing digitized map information, as digitized by 
IOM-TOWELL. 

4. Documentation received from IOM-TOWELL. 

5. Tape containing Landsat classification data. 

6. Documentation of tape contents. 

The requirements for manipulation of these data sets are defined 
In Table 2.2.2 . Each vendor was asked to demonstrate their capability 
In these nine areas, or to indicate how they would meet these require- 
ments If the capability did not exist. In addition to demonstrating 
their systems capabilities, vendors were asked to provide a firm bid 
for installation of the System in Jacksonville. 

During the final System Design Committee meeting in Dallas, Texas 
in early December, vendors capabilities were evaluated. The committee 
was primarily concerned with the vendors capability of meeting the FRIS 
system requirements. Bid information was used by St. Regis staff to 
prepare financial evaluation for St. Regis management. 

2.2.2 Preprocessing Transfer 

LARSYS preprocessing software development is the task which re- 
sulted from a number of FRIS system design meetings beginning in July 
of 1978. As of July 1979, the FRIS system design had progressed to a 
point that the LARSYS preprocessing and analysis software to be trans- 
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ferred to the St. Regis has been determined. LARSYS preprocessing 
software consists of three major processors. The three processors con- 
vert digital Landsat data to LARSYS format, perform systematic geometric 
corrections of Landsat data and register two Images of Landsat data. 

The final task beyond the transfer of software Is the documentation of 
the software. 

Reformatting 

The first preprocessing system of programs converts digital Land- 
sat data to LARSYS format. The functional specifications for this pro- 
cessor required the conversion of input EDIPS "P" Landsat MSS data re- 
ceived In a band Interleaved by line (BIL) format to LARSYS format. 

"P" format refers to EDIPS format computer compatible tape data requested 
as CCT-PM or fully processed MSS data with geometric corrections 
applied and resampled to a map projection. Details of this format may 
be found in the "Manual on Characteristics of Landsat Computer-Com- 
patible Tapes" published by the EROS Data Center in December 1978. 

Preliminary work on the design specifications for the "P" format 
Landsat data processor began in early March 1979. In particular, the 
LARS reformatting group determined that a comprehensive design phase 
would substantially shorten implementation during the programming, de- 
bugging, and documentation phases. Work on the design lasted into 
early July. Every algorithm was defined, program modules specified, 
and nearly all substantive variables and buffer areas were identified 
before programming began. The main routine, along with all subroutines 
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and calling sequences , were thus determined and documented. The program 
destgn Included accomodation for function specifications which would 
support both nearly automatic operations In an operational environment 
as well as multiple options required In a scientific research environ- 
ment. 

Several techniques were utilized to solve this problem. First, 
the basic approach was top-down structured programming. All routines 
have a top to bottom flow of control, and top of calling sequences 
modules were programmed first. As much testing as possible was done 
after completion of each module and assembly of It with previously com- 
pleted modules higher In the calling sequence. The second technique 
was to place a unique or substantive process In a separate module. 
Modules were allocated based on the structured "English" version of 
the processing algorithm (Appendix A-l). The question asked by the 
analyst as he scans this "structured English" program would be what 
processes must occur for this "sentence" or "group of sentences" to 
successfully execute. The answer defined the modules to be programmed. 
Third, the primary programming language was FORTRAN IV utilizing the IBM 
Level 6 or H compilers. Som IBM Assembly Language program includes the 
structured "English" versions of the algorithms used. 

In addition, the Implementation of the entire EDIPS to LARSYS 
programming task was PERT charted. This aided in the management of 
time and resources for the project. Parallel programming efforts 
could then be spotted as well as known or potential bottlenecks. 

Finally, to facilitate the control of the program in the most humanly 


n 


efficient default mode, only three cards are required to execute the 
EDIPS processor. 

Geometric Correction 

The geometric correction processor was the second major system of 
programs. The geometric correction processors original functional 
specifications called for maximum correction of geometric distortions of 
Landsat I data with minimum use of resources. The most important dis- 
tortions thus were corrected. In particular, the data is assumed to 
consist of square 80 meter pixels which are rotated to true north, 
deskewed for the earth's rotation and rescaled for output on a line 
printer with 8 x 10 aspect. In the context of FRIS pre-EDIPS format 
Landsat data may be corrected for geometric distortions. 

In the current FRIS image preprocessing system, this program may be 
utilized to rotate Landsat III data to true north and rescale it if 
necessary. This is especially important considering the number of 
data sources already in true north orientation. Examples are the St. 
Regis Administrative Unit maps. Data in the same orientation is far 
easier to use for the human than data skewed r.- rotated relative to a 
given true north reference data set as the forest AU mentioned pre- 
viously. For example, checkpoints are more readily defined and located 
as part of the image registration process. Relatively minor updates to 
the control card reader to incorporate the rotation-only parameter will 
be required to bring this program to transfer status. Inspection and 
update of program listings, program abstracts, and user documentation 
will also be required. 
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Image Registration 

The last major processor Is the Image registration system. The 
primary purpose of this system Is to register two coincident digital 
images such as two Landsat digital image data sets. The secondary 
purpose is to provide for the registration of any known two-dimensional 
grid to another known or defined two-dlmenstlonal grid. An example Is 
the registration of Landsat data to a US6S quad map. The former has a 
grid X-Y of pixel locations while the latter has a grid of Inches 
horizontally and vertically. Input images are assumed to be in LARSYS 
format. 

The Image registration system software consists of three functional 
sections: 1) the main image registration section, 2) the coincident 

image cross-correlation section, and 3) the multiflt least squares anal- 
ysis. While such software exists to do part of this task, a revised 
processor was desirable to achieve a modern supportable software system. 
The writers of the old system are no longer available and documentation 
for the program Is sparse. The procedure to be followed in this Image 
registration system programming task follows that of the EDIPS pro- 
cessor mentioned above. Once the need for the new processor was 
established, functional specifications were determined. The overall 
goal is to produce a maintainable system which Is modularized, as well 
as documented for program contents, programming techniques and user 
documentation. Furthermore, the latest obtainable registration tech- 
niques for reasonable Implementation. No attempt will be made to dupli- 
cate the Goddard MDP (Master Data Processor). The function of the 
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system, however, Is similar. Two standard registration procedures are 
being utilized to allow more accurate, cost efficient registrations. 

The two Implementation procedures are the use of a cubic poly- 
nomial for the overall registration and the use of blocking with linear 
Interpolation. The first refers to the cubic polynomial whose co- 
efficients will be derived from the MULTIFIT processor. This processor 
uses least squares analysis to derive the best affine, bi -quadratic or 
bl-cubic fits for the checkpoints taken from the respective digital 
Images or known grids as appropriate. With the best equation fit 
determined, normally a bicubic one, the blocking concept will be utilized 
to reduce computation time. 

The concept of blocking during digital Image registration Is a 
moderately complex one. First, the bicubic polynomial for image loca- 
tion Is investigated for rates of change and saddle points by solving 
the first and second derivatives. Utilizing these values one may 
determine the minimum block size within which a bilinear function 
accurately approximates the bicubic one. Block size may be thought of 
as Y lines by Z columns. At least "Z" number of multiple times are 
eliminated from the calculation of each pixel location within the block. 
Only the corner pixel locations of each block need to be calculated in 
full bicubic polynomial mode. The linear Interpolation within the 
block Is relatively fast and predictable with far fewer calculations. 
Should the biquadratic polynomial be the best fit for the data, blocking 
may still be used. However, the reduction in the location calculation 
time will not be as great. In the unlikely event that a linear fit 
will suffice, blocking is not used. 
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Other features of the registration system Include an automated 
cross correlation processor and two forms of pixel gray level inter- 
polation. First the automated cross-correlation processor is an aid 
for acquiring checkpoint locations which are selected from two coinci- 
dent Landsat digital Image data sets. This cross-correlation will be 
accomplished by the Implementation of a numerical integration Image 
correlator. Control of where checkpoints are sought may be by line and 
column Intervals and starting and stopping locations. Alternate con- 
trol may be by a set of arbitrary checkpoints for location after cross- 
correlation. An appropriate Initial transformation will accompany 
either control method. Should this concept not be practical because of 
data dependency problems, manual checkpointing methods will be used. 

The second feature Is a gray level Interpolation method. A gray 
level must be determined for each pixel location In the output grid. 

The nearest neighbor Is the default. The advantage of nearest neighbor 
Interpolation Is that no new data values are created. Classification 
algorithms may use the same statistics before and after registration. 
Cubic Interpolation of pixel gray levels Is the alternative. This 
cubic Interpolation algorithm assumes surrounding pixels input to the 
respective "center" pixel's gray level. The "center" pixel refers to 
the calculated subpixel location outputed from the registration poly- 
nomial. The pixel location Is theoretically subpixel and the level of 
each surrounding pixel to the "center" pixel is determined by which of 
the sixteen subpixel locations Is calculated for the "center" pixel. 

To facilitate the implementation of this third order Lagrange Inter- 
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polation, "center" pixels locations are calculated to one quarter of a 
pixel. Coefficients are pre-supplled In a table for each of the sixteen 
possible "center" pixel locations. The level of calculation Is thus 
^ restricted to simple addition and multiplication. Cubic Interpolation 

of gray levels smooths the visual look of Images. This approach has 
the potential for portraying slightly more accurate subpixel locations 
for given features of the scene. Compared to the nearest neighbor 
i Interpolation technique, the cubic convalatlon approach requires more 

1 computer-resources. 

Preprocessing Documentation 

k 

I 

| Documentation Is the last major effort of the LARSYS preprocessing 

software Implementation. Documentation will consist of three main 

i 

efforts for each of the three processors previously described. The 
three types of documentation will be; 1) program listing documentation, 

f 

2) program module abstracts, and 3) user documentation. The first form 
of documentation has been guided by a standars document produced by the 
reformatting group. These standards expand and clear up details of 
program listing documentation to be followed In the preprocessing 
software. Inputs, outputs, and major variables and arrays are detailed 
[ at the top of each program listing under this standard. In addition, 

; processing procedures are clearly explained through comments In the 

listing. In essence, a new programmer should be able to read the 
comments within the listing and know what algorithm the code Is Imple- 
menting. 
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Programming abstracts are the second form of documentation. 

These follow the LARSYS standard manual. This form of documentation 
normally will be used with the program listing for maximum communication 
to the programmer. 

Finally, user documentation is to be generated. The user document 
should describe what a processor ts used for as well as how to use It. 
Sample control card sets will be Included along with explanations of 
what each set does. It will be coordinated to the extent that a 
minimum amount of Information will be Included regarding where the 
processor Is utilized in the preprocessing of digital images as Landsat 
data. 

Documentation Is a key to the technology transfer of the LARS 
Image processing/analysis system totally known as LARSYS. Good docu- 
mentation although expensive, is necessary to Inform the programmer and 
user. The programs will be more maintainable by less readily knowledge- 
able programming professionals. Over the long term this potentially 
means less total time and expense. To the creator of the documentation, 
the effort means a more thorough knowledge of just what he or she has 
transferred to a fellow programmer in another organization. 

This total LARSYS preprocessing effort has been and continues to be 
a considerable one. Every resource Is being utilized to maintain an 
orderly and effective transfer of the LARSYS preprocessing software In 
a timely fashion. To this point, the schedule o' events has been on 
target. We Intend to keep it that way through its completion. 


2.2.3 LARSYS Transfer 


LARSYS as It currently exists, at Purdue University's Laboratory 
for Applications of Remote Sensing (LARSh consists of 4’ processing 
functions contained In 377 FORTRAN routines and 49 IBM ASSEMBLER 
routines. LARSYS Is not just an Integrated set of computer programs 
designed for the analysis of remote sensing data. It Is an entire 
approach to the conversion of remote sensing data Into Information use- 
ful for monitoring and Inventorying earth resources. Results of the 
Demonstration Phase of FRIS document the utility of the approach to 
Industrial forest management In the southeast. 

The System Transfer Phase of FRIS therefore not only deals with 
the implementation of the LARSYS software, but also the transfer of 
the concept. The FRIS Image processing subsystem will be comprised of 
a subset of the LARSYS ver. 3.1 software package which Is currently 
available through COSMIC. In addition select developmental and experi- 
mental routines, some of which were developed specifically for FRIS, 
will also be transferred. 

As part of this transfer activity, source tapes, program listings, 

m 

users' manuals, system manual , control card references, and program 
abstract* were provided to St. Regis for 23 image processing processors. 
;he following is * brief description of these processors: 

PICTUREPRINT - histograms and display: data In picture form on a 
line printer for each channel selected. 

CLUSTER - using reflectance values from selected channels, groups 
data Into classes and displays the results on a line printer. 
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STATISTICS - calculates transformed divergence between all class 
pairs and performs these calculations for every set of channels 
requested. 

CLASSIFYPOINTS - assigns each pixel In the data to a class, using 
either the maximum likelihood algorithm or minimum distance rule. 
The results are written to tape or disk. 

PRINTRESULTS - using the classification results located on tape or 
disk, prints a map and tabulates the number of pixels c'asslfied 
Into each class. 

IDPRINT - prints most of the Information contained In the MSS data 
header record. 

DUPLICATERUN - duplicates a data run from tape to tape, and 
optionally allows arithmetic expressions to be applied to the data. 
COPYRESULTS - copies classification results from disk or tape to 
another tape. 

LISTRESULTS - prints Information located In the header records of 
the classification results.^ 

PUNCHSTATISTICS - punches a copy of the statistics deck located on 
a classification results tape. 

IINEGRAPH - graphs a line of MSS data on a line printer. 

COLUMNGRAPH - jrajihs a column of MSS data on a line printer. 
HISTOGRAM - histograms data and produces a deck of the histogram 
Information. 

GRAPHHISTOGRAM - on a line printer, displays the histogram produced 
by PICTUREPRINT or HISTOGRAM processors. 
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ECHO - extracts and classifies homogeneous objects as If they were 
single pixels. 

MERGESTATISTICS - combines more than one statistics deck Into a 
single deck. 

BROWSt - searches the LARSYS runtable and prints the list of MSS 
data runs which satisfy the us—'s request. 

RATIO - using the mean vectors of classes In a statistics deck, 
calculates and prints the ratio of the values for the specified 
channels and the sum for each class. 

GOATA - histograms and displays data In picture form on a 
printer/plotter for each channel selected. 

GRESULTS - using the classification results, prints a map on a 
printer/ pi otter and tabulates the results. 

BIPLOT - produces a blspectral plot of classes contained In a 
statistics deck. 

CHANGE - compares two classifications results over the same area 
for purposes of change detection. 

The above processors represent 42,000 lines of FORTRAN, 5500 lines of 
Assembler, a.nd 1500 lines of CMS EXEC language. The conversion effort 
for this software consists of 1) converting from FORTRAN G to FORTRAN 
H, 2) removing operating system dependencies from the Assembler, and 
3) convertin’! ths CMS EXEC language into JCL. During this past quarter, 
the first two or these were accomplished, although the software will 
not be able to be tested until the third Is completed. 


A - • 
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2.2.4 Programming Additions 

During this quarter, the fnttlal coding and debugging of a new 
L.ARSYS processor was completed. TMs processor Implements a polling 
algorithm which scans a classification results file a cell (or group 
of points) at a time, and assigns the majority class number from that 
cell to all points in the cell. The purpose of this Is to eliminate 
isolated points from the classification map in order to provide an out- 
put product that matches the naming conventions of the user. In this 
case the St. Regis Paper Company. This processor is currently being 
tested and enhanced. Documentation and sample maps will be included 
in the next quarterly report. 

2.3 Applications Test 

In mid- July 1979 the FRIS APT Manager requested that the project 
staff undertake an operational test of the remote sensing technology. 
Specifically, the project staff was to provide a change classification 
of a tract of land in Baker County, Florida (figure 2.3.1) that St. 
Regis had recently acquired. 

The now acquisition, the Knabb Tract, encompasses approximately 
40,000 acres of land and is ecologically similar to the Fargo test 
site. St. Regis staff were of the opinion that the timber removals 
had been extensive in recent years. Furthermore, they felt that re- 
movals were especially extensiye since 1977. 

The application test was designed to address the feasibility of 
using Landsat classified data to: 



cation of Baker County, Florida and the Knabb Tract 
•plication Test Site. 
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1) Evaluate the areal extent of the standing timber resource, from 
1979 data, and 

21 determine the change in standing timber that was detected by 
Landsat that occurred between 1977 and 1979. 

) Our goal was to meet these objectives and to provide classification 

results hy 1 November 1979. Timing was an important criteria to this 

i application test because if the data could not be: 

I o Acquired 

l 

■ o Preprocessed 

o Classified, and 

i 

o Final products available 

by the deadline, than the timeliness of the technology would be 
seriously questioned. In November the window for photographic data 
collection opens, and this tract was flown. If Landsat information was 
not available by the time the photography was collected and interpreted 

1 than the utility of Landsat as an industrial forest management tool 

will be seriously questioned. 

Two critical elements to the successful completion of this appli- 
cations test were the data preprocessing and change classification 
tasks. Descriptions of these activities follow. 

Knabb Tract Data Preprocessing 

The primary preprocessing task involved the registration of two 
Landsat frames to a 1;24,00Q scale base with property boundary infor- 
mation merged with the Landsat imagery. Although one of the Landsat 
frames (21050-14515) was available in-house, a portion of the frame 
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necessary for the preprocessing was destroyed during an earlier process, 
necessitating the reordering of the data. The second data set, Landsat 
frame 21482*15101 was not expected to arrive until 28 September 1979. 
Both Landsat frames would have to be reformatted to LARSYS Version 3.0 
format and geometrically corrected (systematic removal of first order 
distortions) to a scale of 1:24,000 with a line printer aspect 
(10X:8Y). 

As one of the intents of the test was to determine the timeliness 
of Landsat in providing land cover information, it was important to 
complete the preprocessing activity in minimal time. Using PERT 
planning, a probable completion date was estimated as 12 October 1979. 
This date was based upon a starting date of 1 August 1979 and receipt 
of the February 1979 Landsat data on 28 September 1979. Actual com- 
pletion of the preprocessing task was 11 October 1979. 

Map digitization was performed by St. Regis Southern Timberlands 
Division personnel. The original intent was to digitize a 1 inch to 
1 mile map in the hope of eliminating the editing problem of digitally 
reconnecting the maps. However, the boundaries to be digitized were 
drawn on 1:24,000 scale US6S quadrangle maps. A determination was 
made that accuracy would be lost by transferring the boundaries to the 
1 inch to 1 mile map and then rescaling the data back to 1:24,000. A 
decision was made to digitize the boundaries directly from the USGS 
1:24,000 scale maps and join the maps together digitally. This final 
method worked very well with no unanticipated problems. A portion of 
the digitized data Is shown in Figure 2.3.2 . 
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At the same time the maps were being digitized and the digital 
boundary information edited the 7 December 1977 data was reformatted 
to LARSYS format and geometrically corrected. After completing the 
digitizing, 14 checkpoints were located hetween the 7 December 1977 
data and the USGS quad maps using the LARSYS IMAGEDISPLAY program. 

The 14 control points were run through an affine (6 parameter non- 
conformal ) least squares fit. The resulting transformation function 
exhibited a line error of 0.708 (rms) and a column error of 1.032 (rms). 
The following first order distortions were corrected by the transforma- 
tion of the systematically corrected 7 December 1977 data to the 
1:24,000 USGS map coordinates: 

Scale X 1.0152 
Scale Y 1.0000 
Rotation 0.326 degrees 
Skew 0.0299 degrees 

At this point, both the digitized map boundaries and the Landsat data 
were in the same reference coordinate system. 

The next preprocessing step was to actually create the ownership 
information in grid form. This was accomplished by "rasterizing" the 
vectored digital boundary data. Some editing of an intermediate file 
is normally required when the boundaries to be rasterized are of 
regular rectangular polygons. This was the case for the Knabb tract 
although minimal editing of the intermediate file was required. The 
final result was a tape in LARSYS format containing the precision (map) 
registered data with an auxiliary data channel containing ownership 


26 


information. All data outside of the ownership was set to a null 
value (hex QO). 

The second Landsat data set (21482-15101, 12 Fehruary 1979} 
arrived on September 28. The tape received was created in the newer 
EDIPS format rather than in the expected X format. The tape was re- 
formatted using the new EDIPS to LARSYS reformatting software. 

As the second data set was already corrected by NASA to eliminate 
radiometric and geometric distortions, it was necessary to correct only 
for rotation, scale, and aspect. The EDIPS corrected tapes are regis- 
tered to a Mercator projection (either Hotine or Space Oblique} using 
ground control. The resulting scale of this data set is approximately 
1 :179.52.7 with each pixel representing a ground resolution of 57m 
square. As the current geometric correction processor is designed to 
correct for pixel size and skew, corrections already applied to the data 
by the NASA process, a transformation to correct for rotation, image 
scale, and pixel aspect using the image registration system was developed. 
The method for performing this type of correction through the image 
registration system is described in Appendix A2. 

The final step was to register the corrected second Landsat scene 
to the December 1977 data. A total of 185 checkpoints were located 
between the two images using the numerical autoucrrelator of the image 
registration system. An average correlation coefficient of 0.69 was 
obtained through 270 correlation attempts between the second channel of 
each scene. The average error between the predicted and observed 
checkpoint location was 0.67 pixels. The checkpoint pairs were then 
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run through a biquadratic least squares fit. All control points were 
accepted with rms errors of 0.099 in the line direction and 0.283 In 
the column direction. The following first order distortions were 
corrected by registering the corrected EDI PS Landsat scene to the map 
reference grid: 


Scale X 

0.9999183 

Scale Y 

1 .0006683 

Rotation 

-0.1515851 degrees 

Skew 

0.075 degrees 


The Knabb Classification 

The acquisition by St. Regis of the Knabb tract provided an oppor- 
tunity to extend the classification procedures into an unknown area - 
one for which no photography or forest cover type information was 
available to the analyst to aid in defining a classification training 
set. The Knabb property is located about 50 miles southeast of the 
Fargo, Georgia. The Fargo test site had been classified in the first 

year of the project. Both areas are included in the same Landsat data 

faames and the forest cover types were known to be similar at both sites. 
Since there was interest in finding how quickly such a classification 
could be produced, the initial plan was to use the existing training set 
from the Fargo area to classify the* Knabb site. 

The most recent data set from these sites had been collected on 
December 7, 1977, although a data set collected on December 30, 1976 

was also available. After examining the data it was found that in the 

1976 data the Knabb site was about half obscured with clouds and most of 
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the remaining half was In cloud shadow. The 1977 data was cloud- free 
and otherwise of excellent quality. However, no single data training 
set was available for the 1977 data since It had been used only In a 
bitemporal analysis. 

In order to save time the December 7, 1977 data was classified 
with the December 30, 1976 training set. Normally this difference In 
data sets would have caused serious. If not Insurmountable, data cali- 
bration problems. In this case, however, the dates of data collection 
were both In December and were very near to the data of minimum sun angle. 
This, plus the fact that the weather condition were ideal over the 
training area In 1976 and the Knabb tract in 1977 allowed the use of the 
1976 training set with 1977 data without calibration. The only signifi- 
cant classification problem was found to be the lack of a training 
class to represent the clear cut/site prepared areas which were present 
at the Knabb site but did not occur at Fargo. This deficiency was 
quickly corrected by adding two training classes generated from the 
Knabb 1977 data to represent these cover type conditions. A summary Is 
shown In Table 2.3.1 . 


Table 2.3.1 Area statistics for the Knabb Tract calculated from a 
classification of December 7, 1977 Landsat data. 


Cover type 

Acres 

Hectares 

Percent 

Pine 

22,723 

9,200 

52.2 

Pine/Hardwood 

10,916 

4,420 

25.0 

Slash/Cypress 

8,275 

3,350 

19.0 

Nonstocked 

1,521 

616 

3.5 

Wet lands 

122 

49 

0.3 


43,557 

17,635 

100.0 
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After this Initial classification was completed a new data set was 
received. This data set, collected on February 12, 1979, was overlain 
onto the 1979 data. Property boundary lines were digitized and added 
to this data set. A separate analysis was carried out using the previous 
classification augmented with Information gathered during the field 
checking as training aids. The data quality was not nearly as good as 
that of the two previous sets. The 0.6-0. 7 mlcormeter band was unusable 
due to severe banding. The classification was done with the three re- 
maining bands. A summary of the classifications is shown in Table 2.3.2 . 


Table 2.3.2 Area statistics for the Knabb Tract calculated from a 
classification of February 12, 1979 Landsat data. 


Cover type 

Acres 

Hectares 

Percent 

Pine 

25,487 

10,319 

58.5 

Pine/Hardwood 

5,972 

2,418 

13.7 

Slash/Cypress 

10,959 

4,437 

25.2 

Nonstocked 

719 

291 

1.7 

Met lands 

420 

170 

0.9 


43,557 

17,635 

100.0 


In the two-year interval between the two data collections, several 
areas were cut and planted or were being prepared for planting. The two 
classifications were compared with the CHANGE processor to find and 
Identify these areas and those results are shown in Table 2.3.3 . The 
two original classifications are shown In figures 2.3.3 and 2.3.4 and the 
change map is shown in figure 2.3.5 . 
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Table 2.3.3 Area statistics for the Knabb Tract showing changes In 

ground cover which occurred between the two classifications. 


CRange 

Acres 

hectares 

Percent 

New Plantation 

7,686 

3,112 

17.7 

Harvested 

484 

196 

1.1 

No change 

30,634 

12,403 

70.3 

Unidentified 

change 

4,753 

1,924 

10.9 


43,557 

17,635 

100.0 


2.4 Management 

The management activity oversees day to day project operations and 
Is responsible for all technical and fiscal project reports. Status of 
all major System Transfer Tasks are shown In Exhibits 1 through 4 of 
Appendix B. 

In addition to the operational project functions, special support 
studies, like the Knabb Applications Test (reported In section 2.3) are 
monitored from the task. One project that Is especially noteworthy in- 
volves the application of regression analysis to Landsat MSS data to 
predict age of managed southern pine plantations. More precisely the 
ratio of the Infrared to visible band responses are assumed to be 
affected by stand occupancy, which is reflected In crown closure. As 
stands mature, individual tree crowns occupy a greater proportion of the 
site (figure 2.4.1). The increasing crown closure affects the ratio, 
which In preliminary tests corresponds well to a measure of age. 
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Figure 2.4.1 


This Is a conceptual representation of the biological 
growth response over ttae. The b.stc hypothesis of . r.tlo 
evaluation assumes that as a plantation matures the reflec- 
tance of the tree crowns will be the dominant factor affect* 
Ing the calculated ratio. Therefore, this measure of crown 
closure can be related to stand age. 


i 
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Knabb and Picayune Ratio Results 

The ratio of IR che nne Is to visible channels from December 1977 
Landsat data for the Knabb and Picayune tracts was used to predict the 
age of selected pine fields. The exact ratio used, the method of 
picking pine fields, the analysis used to predict the fields' ages and 
the results of these predictions are outlined below. 

The exact data ratio generated was as follows: 

ratio - 40.0(C3 ♦ C4)/(C1 + C2 ♦ 0.1) 

where 

Cl » channel 1 
C2 ■ channel 2 
C3 ■ channel 3 
C4 ■ channel 4 

The multiplier 40.0 and the constant 0.1 were needed to enhance the 
range of and Information in the data, and to prevent a divisor of zero. 

The Knabb tract was first categorized Into pine and non-pine 
classes. From this classf^lcatlon fourteen fields of seemingly homoge- 
neous pine were selected and the average ratio for each field was 
determined. Due to the proximity of this tract to the Fargo test site 
and their similar physiography, a regression equation developed for 
Fargo was used to predict the ages of the selected Knabb fields. Four 
of the Knabb fields were dropped from further analysis. Two of these 
discarded fields were accidently picked outside the Knabb boundaries 
and the other two dropped fields wer? Inaccessible for checking ground 
truth. Of the ten pine fields left, a ground Inspection of tie area 
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established that (1) all ten fields were pine, and (2) nine of the ten 
fields had ages within the ninety percent confidence Interval for each 
predicted age. Ages were derived by taking Increment cores and 
counting growth rings of randomly selected dominant trees. 

Table 2.4.1 shows the preliminary results obtained by applying the 
Fargo prediction equation to the ratio calculated from landsat data over 
the Knabb Tract. 


Table 2.4.1 Preliminary results for ten pine plantations in the Knabb 
Tract. 


rge Measured 



Prediction Equation: 

1og 10 (AGE ♦ 1) - -9.333088 + 5.567559 log 1Q (ratio) 

Similarly, ten pine fields were chosen from the Picayune Test Site 
and the predicted age of each field was determined from *ii equation 
developed specifically for the Picayune data. The age ~ ■ each field 
was verified by ground investigation. Table 2.4.2 shows the preliminary 
resul ts. 
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Table 2.4.2 Preliminary trials of an age prediction equation using 
Landsat ratio values for Picayune, Mississippi. 


Field 

1 

2 

3 

*4 

5 

6 

7 

8 

*9 

10 


Ground 

Verified Age 

Age Pred 
from Lands 

15 

17 

18 

9 

14 

10 

26 

6 

2 

3 

14 

10 

13 

18 

0 

1 

26 

9 

14 

16 


Prediction Equation: 

log 10 (AGE + 1 ) = -4.913634 + 3.218647 log ]0 (ratio) 

♦both these fields had actual ages beyond the range of the regression 
equation. Of the ten Picayune fields checked, two (fields 4 and 9) 
fell outside the 90% confidence interval for the predicted age. 

Another application of the generated ratio channel was a classi- 
fication of the Knabb area done solely with the ratio channel 
(LEVELCLASSIFY) . Analysis done on the Fargo test site revealed the 
fact that the average ratio of hardwood fields in winter data fell below 
the average ratio of pine fields. Hence using the ratio intervals 
developed on the Fargo test site, the Knabb tract was classified into 
hardwood, young pine (less than 15 years old), and old pine (15 years old 
or over). 

Since the levels for the level classifier were determined using 
averages over fields, these levels did not apply directly to classifying 
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Figure 2.3.3 Classification of December 7, 1977 Landsat data for the Knabb 
Tract. Area statistics for this classification appear in 
Table 2.3.1 . 
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pixels. Also the levels were determined on another site causing even 
more inherent error in this classification. The ten pine fields used 
in Table 2.4.1 and four hardwood fields were used to test the accuracy 
of this classification. The results are presented in Table 2.4.2 . 


Table 2.4.3 Classification performance for a LEVF.LCLASSIFY approach 
using Landsat ratio input for the Knabb Tract. 

Classes 

Classification accuracy 

young pine 

45.8 

old pine 

57.4 

hardwood 

60.0 

overall accuracy = 57.4 



Hence the ratio of IR to visible Landsat channels has shown useful- 
ness in predicting the ages of pine fields even over areas with no 
ground truth. 


Preliminary results using the generated ratio as a classification 
channel, however, has shown questionable usefulness. This does not 
preclude further investigation of levels classification technology. 

The levels classifier is significantly faster than a maximum likelihood 
per point approach and could therefore be beneficial for "first look" 
evaluations of large areas. Additional investigation into the appli- 
cation of this approach will be pursued. 


APPENDIX A-l 


Example of Structured English Module 
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Structured "English" Program 


PROGRAM lndsup 

(♦ Landsa t to LARSYS program •) 


DECLARATIONS 

control. cards : file; 

end-of-filo s file-condition; 

error-level, abort : error. level-indicator; 
start, stop : tine? 

dummy, clock, totcpu, vircpu : computer. usage? 


BEGIN 

CALL cptinc (dummy , clock, totcpu, vircpu); 

CALL getinc (start) ; 

REPEAT 

CALL lrdint; ( # initialize *) 

CALL lndrdr ; ( # read control card file *) 

IF NOT (cnd-of- file ON control-cards) THEN 
IF error-level < abort WEN 

CALL lndctl ; (* rofcrrcvat Landsat *) 

CALL InJwup (* wrap-up reformatting # ) 

END IF; 

CALL it.dsum ( # summarize the job •) 

ENDIF 

UNTIL end-of-filu ON control. card3 OR 
error-level >=* abort; 

CALL usage (clock, totcpu, vircpu); ( # print corrputer time used •) 
CALL getimo (stop) ; 

CALL lndf 17 (* print FORM- 17 *) 


END 
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APPENDIX A-2 


Instructions to Correct EDIPS Corrected Data Using the LARS Image 


Registration System 
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Instructions to Convsrt Linear Parameters to 4-Parameter Non-Conforral 
Transformation for use by the LARS Image Registration System. 

(This algorithm is for use with Landsat-EDIPS formatted geometrically 
corrected tapes). 

1. Set up the distortion matrix M such that 


where 


M - M 1 M 2 M :) M 4 M : 


*1 


1.0 o.ol 


0.0 1.0 


for pixel scale (assumes 57m pixel) 


*2 


cos 0 sin 0 

'Sin 0 cos 0 


for rotation, where 0 is the rotation 
required to a north-south orientation. 


1.0 o.o' 
1.0 0.0 


0 ■ sin 


n 


(■ 


cos ( 80.985) m 
cos $ 


where $ is run center latitude 
for skew (already corrected by 


EDIPS) 


M. 


0.8 0.0 

0.0 1.0 


for a line printer aspect of 8 lines 
per inch, 10 columns per inch 


# 


*5 


1.3368421 O.o] 
0.0 1.3368421 


for scale correction to 1:24000 from a 
scale of 1:17952.7 


■sr 




9 


i 

1 
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M «* 




1.0694737 cos Q 1.336842) sin 0 


-1.0694737 sin 0 1.3368421 cos ©j 


This distortion matrix is set up for the transformation 
V - M X 

where X is the original coordinates and 
7 is the new coordinate matrix. 


2. To put this into a form suitable for the registration equation: 

Y • X ♦ A 

we must use the form 

Y » A X 

Expanding: Y-AX“X. + A 

or Y^ * a n x i + a i2 X 2 * ax l + bx 2 + c + x i 

*2 ‘ *21 X 1 * ®22 X 2 ' d *l * "2 * f * X 2 
let c ■ f ■ 0, 

*n«l * * l 2 x 2 * <* * 11 X 1 * bx 2 


*21 X 1 * *22 X 2 - dx l * <* * llX i 


or. 


*u • * * 1 
*12 * b 
*21 

*22 ' ' * 1 


or, 


a + 1 b 
d e + 1 
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3. To convert A in terms of overlay line (CLD) and column (CJD) 
coefficients of the image registration system, it is only 
..acessary to observe 

a * CJD 3 = 1.0694737 cos 0 - 1 
b = CJD2 = 1.3368421 sin 
d « CLD 3 =-1.0694737 sin 
e = CLD2 = 1.3368421 cos 0 - 1 
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I MAGE RE GISTRA TION PE RFOI ^'JC^ REPORT 


DATE: 9/25/79 


REFORMATTER: Smith 


W P- 556-PR 
MACT- 796 


CPU- 

INPUT 


RUN A 77010200 LINES 1,600 COLS 75,435 CHAN _ 

RUN B 77010201 LINES 1,750 COLS500, 1100CHAN All 

BUFFER USED: 


OUTPUT 


RUN 77010202 TAPE 2668 FILE 2 

INTERPOLATION: NN 

DATA TRANSFORMATION INFORMATION 


TRANSFORMATION 
LINE RMS 
COL RMS 
N 

DISTRIBUTION R 


0.70806 

1.03212 

14 

0.99 


FORWARD 


^backward) 


ORDER 1 


Registered to USGS Quad Maps 


CORRELATION NA 


RATE OF ACCEPTANCE 


ACCEPT REJECT RATE 


AVERAGE RHO 

AVERAGE EUCLIDEAN ERROR 
TRANSFORMATION COEFFICIENTS 


LINE COEF. 1: 

32.494452 


COL COEF. 

1: 

543.055448 

2: 

0.015229 



2: 

-0.005780 

3: 

0.005092 



3: 

-0.013307 

4; 




4: 


5: 




5: 


6: 




6: 


THE FOLLOWING 

FIRST ORDER 

DISTORTIONS 

WERE CORRECTED: 



SCALE 

X: 1.0152 

Y: 

1 

.0000 


ROTATION 

0.326 

(DEGREES) 



SKEW: 

0.0299 

(DEGREES) 



IMAGE REGISTRATION PERFORMANCE REPORT 
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DATE: 10/5/79 


REFORMATTER: Smith 


WP« 572-GC 


MACT* 796 


INPUT 


CPU* 


RUN A 79000300 LINES 762,1500COLS 1200,1698 CHAN none 

RUN B 79000200 LINES 600,1700COLS 1000,2150 CHAN 1-4 

BUFFER USED: 1,500,000 

OUTPUT 

RUN 79000201 TAPE 2836 FILE 1 

INTERPOLATION: 


DATA TRANSFORMATION INFORMATION 


TRANSFORMATION 
LINE RMS 
COL RMS 
N 

DISTRIBUTION R 


FORWARD < r^ACKWAR5) ORDER 1 


NA: systematic distortions corrected for rotation, 

scale , aspect only 


CORRELATION 


NA 


RATE OF ACCEPTANCE ACCEPT REJECT RATE 

AVERAGE RHO : 


AVERAGE EUCLIDEAN ERROR : NA 

TRANSFORMATION COEFFICIENTS: 


LINE COEF. 1: 

0.0 

COL COEF. 1: 

0.0 

2: 

0.314711 

2: 

0.2422435 

3: 

-0.1937948 

3: 

0.0517688 

4: 

5: 

6: 

THE FOLLOWING 

FIRST ORDER DISTORTIONS 

4: 

5: 

6: 

; WERE CORRECTED: 



1.0694737 

SCALE X: 1-.2400C LP Y: 1:24000 LP 

ROTATION 10°. 44 (DEGREES) 

0°,0 (DEGREES) 


SKEW: 


IMAGE REGISTRATION PERFORMANCE REPORT 
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wp- 583- IR 

DATE: 10/11/79 REFORMATTER: Smith 

MACT* 796 


CPU* 

INPUT 


RUN A 77010202 LINES 1,600 COLS 1,428 CHAN 1-4 

RIM B 79000201 LINES 1,739 COLS 1,504 CHAN 1-4 

BUFFER USED: 7246 bytes 

OUTPUT 

RUN 79000202 TAPE 2848 FILE 1 

INTERPOLATION : NN 

DATA TRANSFORMATION INFORMATION 


TRANSFORMATION 
LINE RMS 
COL RMS 
N 

DISTRIBUTION R 


FORWARD 



0.099 

0.283 

185 

NA 


distribution appears excellent 


ORDER 2 


CORRELATION 


RATE OF ACCEPTANCE 


ACCEPT 185 REJECT 85 RATE 


AVERAGE RHO 

AVERAGE EUCLIDEAN ERROR 
TRANSFORMATION COEFFICIENTS 


0.69 

0.67 


LINE COEF. 


129.20298180 
-0.00008529 
-0.00133339 
-0.00000039 
0.00000173 
0.00000115 

THE FOLLOWING FIRST ORDER DISTORTIONS WERE CORRECTED: 
SCALE X: 0.9999183 Y: 

ROTATION -C. 1515851 (DEGREES) 

SKEW: 0.075 (DEGREES) 


COL COEF. 1 

2 

3 

4 

5 

6 


66.39532730 

0.00264544 

0.00066831 

-0.00000293 

-0.00000105 

0.00000014 

1.0006683 


A -4' ;'« f U! I JK iKlf A 



'•*!!(. it liiCLlOlAH fck 
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Appendix B 

The Exhibits that appear in this Appendix are updated versions of 
the milestone and timeline charts that were published in the first 
Phase III quarterly report. The Exhibit numbers to the following 
Phase III task activities: 

1. Technology Transfer 

2. LARSYS Transfer 

3. Preprocessing Transfer 

4. Management 


54 


Exhibit 1 

FRIS III Timeline Chart 



Calendar Year J 

Task: Technology Transfer 

1979 ! 

1980 i 

m 

mm 

10/1 - 
12/3 

1/1 - 
3/31 

4/1 - 
6/30 

Activity: 






A, Training 

1. Short Courses 

2. Workshops 

3. fkno- I nterpretation 

Snort Course 

▼ 

▼ 



mum m ■ • < 

,.S7 

B. Consultation 

m . 


n 

.~y 


C, Ddcufentation 

1. LARS User Documentation 

2. NCC User Documentation 

D. Terminal Operations 

i 


▼ 


S7 

W - 




y 

! 



1 

i 





▼ 

tL 


planned start of activity 
actual start of activity 
duration of activity 
progress toward activity completion 
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Exhibit 2 

FRIS III Timeline Chart 


Calendar Year 


Task: LARSY5 TRANSFER 


Activity: 

A. Planning 

B. Transfer 

C. Consultation & Debugging 

D. DoCUTtNTATION 

E. Test & Evaluation 



planned start of activity 
actual start of activity 
duration of activity 
progress toward activity completion 
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Exhibit 3 


FRIS III Timeline Chart 


Task: PREPROCESSING TRANSFER 


Activity: 

A. ftjWNING 


B. Program Refinement 

1. Landsat 3 Reformatting 

2. Geometric Correction 

3. Image Registration 


C, Program Transfer 

D. Consultation & Debugging 

E, DOOIENTATION 

F. Test & Evaluation 


G. Sipport Activities 

1. Landsat 3 Data Evaluation 

2. FRIS Map Coordinates Definition 

3. ^formatting Operations 

Procedures 


I Calendar Year 

1979 ' I 

1980 1 

4/1 - 
6/30 

7/1 - 
9/30 

10/1 - 
12/3 

1/1 - 
3/31 

4/1 - 
6/30 

mr. 


^ 

f S7 

.-Sr 

* 

• 

w 


W .. 



V._ 

i 






'C 



X7 


V 




_J7 

_S7 











r 

▼ 


planned start of activity 
actual start of activity 
duration of activity 
progress toward activity completion 


OM 
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Exhibit 


FRJS III Timeline Chart 


Task: MANAGEfENT 


Activity: 

A. Reporting 

1. Informal Monthly Status 

2. Monthly Fiscal 

3. Quarterly Progress 
A. Semi-Annual fe views 


Calendar Year 



7/1 - 
9/30 



1/1 - 4/1 - 

3/31 6/30 


B. Information Dissemination 

C. Cost Evaluation 

D. Special Projects 

1, Classification Accuracy 

Evaluation 

2, Ratio Evaluations 

3, Knabb Application Test 


planned start of activity 
actual start of activity 
duration of activity 
progress toward activity completion 








